Realizing and Detecting the Haldane's Quantum Hall effect with Ultracold Atoms 
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We design an ingenious sciieme to realize the Haldane's quantum Hall model without Landau levels by using 
ultracold atoms trapped in an optical lattice. Three standing-wave laser beams are used to construct a wanted 
honeycomb lattice, where different on-site energies in two sublattices required in the model can be implemented 
through tuning the phase of one laser beam. The staggered magnetic field is generated from the light-induced 
Berry phase. Moreover, we establish a relationship between the Hall conductivity and the atomic density, 
enabling us to detect the Chern number with the typical density-profile-measurement technique. 
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The quantum Hall effect (QHE) in two-dimensional 
electron systems is one of the most pecuUar quantum- 
mechanical phenomena observed in nature. The QHE is usu- 
ally associated with a uniform external magnetic field, which 
splits the electron energy spectrum into discrete Landau lev- 
els (LLs). When the Fermi energy lies in the gap between two 
LLs, the Hall conductivity in units e^/h is accurately quan- 
tized to an integer. The precise quantization of the Hall con- 
ductivity was explained by Laughlin based upon a gauge 
invariance argument, which is fundamental to the picture of 
edge states proposed by Halperin On the other hand, 
Thouless, Kohmoto, Nightingale and Nijs (TKNN) inter- 
preted the Hall conductivity as the topological Chern number 
of the U{1) bundle over the magnetic Brillouin zone of the 
bulk states. 

Twenty years ago, Haldane showed in principle that a QHE 
may also result from breaking of time-reversal symmetry 
without any net magnetic flux through a unit cell of a periodic 
two-dimensional (2D) system, where the electron states retain 
their usual Bloch state character jst). In his work, Haldane 
constructed a tight-binding model on a honeycomb lattice in- 
cluding a complex second nearest-neighbor hopping integral. 
The honeycomb lattice consists of two triangular sublattices A 
and B with different on-site energies AI and —M, as shown 
in Fig. la. For M ^ 0, the inversion symmetry is broken 
and the lattice possesses the point group C^v symmetry. A 
periodic vector potential A(r) is applied to the lattice, given 
that the total magnetic flux through each unit cell vanishes, 
i.e., the first-neighbor hopping integral t is unaffected. The 
second-neighbor hopping integral t acquires a Peierls phase 
factor cxp(ie J A ■ dr/fi), where the integration is along the 
hopping path. The Hamiltonian of the model is written as 

H = ^{ta]bj + H.c.) + M{a]aj - b]bj) 

+ ^ /e*'^^'(a|a, (1) 

where and hi are the annihilation operators on site Ri in 
sublattices A and B, respectively, i^ji is the accumulated 



Peierls phase from site j to its second neighbor /, which is 
assumed to take the form ipji — ±Lp. The hopping direc- 
tions for which Lpji = +tp are shown in Fig. la. The most 
interesting and unique feature of the model lies in that the 
phase of the system can be changed from a normal insula- 
tor to a Chern insulator by the simulation of parity anomaly 
il S S S H ES El- However, it is extremely hard to 
realize the Haldane's model experimentally in ordinary con- 
densed matter systems because of the unusual staggered mag- 
netic flux assumed in the model. 

On the other hand, the technology of ultracold atoms in an 
optical lattice provides a perspective approach to explore rich 
fundamental phenomena of condensed matter physics lll2lll3 . 
1411 . In particular, how to realize theOHE with cold atoms has 



attracted considerable interest ! 15l Il6l |17L [18l [19|]. Neverthe 



less, the atomic QHE has not been observed yet, mainly due to 
challenges in both the realization and detection of the atomic 
Hall effects. Although an effective magnetic field for neutral 
atoms can be simulated either by rotating the atoms ll20ll or by 
laser-induced Berry phases ll2U l22l l23l 12411 . the strong mag- 
netic field region required for QHE has not been reached yet 
in experiments. For the rotating method, the system is close to 
the point at which the centrifugal potential cancels the exter- 
nal harmonic trap, and the atoms may fly apart at the rotation 
speed required by QHE 1 3 13] ■ For the laser-induced Berry 
phase approach, the cold atoms moving in a spatially varying 
laser field feel an effective gauge potential! ill 
but the region of the strong uniform field is rather small for 
two typical counterpropagating Gaussian laser beams. In ad- 
dition, the detection method for cold atoms is very differ- 
ent from that for condensed matter systems; especially the 
widely used technique for QHE based on the transport mea- 
surements is not workable for atomic QHE. In this Letter, we 
design an ingenious scheme to realize the Haldane's quan- 
tum Hall model without LLs by using ultracold atoms trapped 
in an optical lattice. We work out a distinct method to con- 
struct the honeycomb lattices that have different on-site ener- 
gies by three standing-wave laser beams. Although it is still 
hard to achieve a strong homogenous magnetic field required 
by the conventional QHE in atomic system, we may evade 
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FIG. 1; (Color online) (a) Illustration of the Honeycomb lattice struc- 
ture of graphene, where open and solid circles represent sites in sub- 
lattices A and B. ai — (ia, — -^a), a2 = (^a, -^a) are the 
unit vectors of the underlying triangular sublattice. si, S2 and S3 
are three vectors pointing from a B site to its three nearest-neighbor 
sites, (b) and (c) show the contours of the potential V for x ~ 2n/3 
and X = 397r/60, respectively. The vertical (horizontal) axis repre- 
sents ykg /tv {xk^ /tt). (d) Contours of the magnetic field defined by 
Eq.lH. 



this bottleneck since LLs are not necessary in this unconven- 
tional QHE. We elaborate that the staggered magnetic field, 
which is hard to generate in condensed matter systems, may 
be rather easy to set up by other three standing-wave laser 
beams. In this scenario, different on-site energies in two sub- 
lattices can be easily adjusted through tuning the phase of one 
of the laser beams, and thus the whole phase diagram|5] in- 
cluding the exotic topological phase transition predicted by 
Haldane may be revealed experimentally. Furthermore, based 
on (2 + 1) -dimensional relativistic quantum mechanics cal- 
culations, we establish a direct relationship between the Hall 
conductivity and the equilibrium atomic density, such that the 
famous topological Chern number may be experimentally de- 
tected with the standard density profile measurement used in 
atomic systems ||l2i. 

Let us first consider single component fermionic atoms 
(e.g., ""^K, ''Li, etc.) in a 2D honeycomb lattice HQ, which 
can be realized by three detuned laser beams. A detuned 
standing-wave laser beam will create a potential in the form 
Vq sin^(kQ • r), where Vo is the potential amplitude and kg 
is the wave vector of the laser. To generate the honeycomb 
lattice with different on-site energies in sublattices A and B, 
the three laser beams with the same wave length but different 
polarizations are applied along three different directions: 
and ^Cj. ± iey, respectively. The potential is thus given by 
V = Vb[sin2(a+ + f ) + sin2(yfc^ + f ) + sin^ia^ - f )], 
where a± = y/ixk^ /2 ± yfcg /2. The potential contours are 
plotted in Fig. lb and Fig. Ic. Under these conditions, atoms 



are trapped at the minima of the potential, forming a honey- 
comb lattice. An amazing feature here is that the different site- 
energies of sublattices A and B is controllable by the phase of 
laser beam x- For instance, we get exactly the honeycomb 
lattice with the same on-site energies {M = 0) for x = |7r, 
as shown in Fig. lb, while the two sublattices have different 
on-site energies { AI ^ 0) for x 7^ 

Now we elaborate how to simulate the staggered mag- 
netic field in the Haldane's model. Since the net flux per 
unit cell vanishes, the vector potential applied to the lat- 
tice must be periodic. Such magnetic fields can be cre- 
ated by Berry phase induced from two opposite-travelling 
standing-wave laser beams ll22ll . For the two laser beams 
with Rabbi frequencies 51i = fio sin(yfc^ + £)e'^'^i and 

^2 = r2ocos(yfc2 + |-)e^'^'^f, the effective gauge poten- 
tial is generated as Ai(r) = sin(2yfc|')e2; 12211 . Here, 
fcf = cos 9 and fcl" = sin 9 with k^ the wave vector of 
the laser and 9 the angle between the wave vector and the 
axis. We emphasize that the choice of wave vector fcj of the 
laser beams must be a multiple of in order to be com- 
mensurate with the optical lattice. We take fcl" — ^^f^- The 
Peierls phases for the nearest-neighbor hopping in Fig. la are 
V12 = ^61 = = -<^45 = Va and 1^23 = y^se = 0. For 

the next-nearest-neighbor hopping integrals, which are inte- 
grated on a period of the vector potential, the coiTesponding 
accumulated phases are cpi^ = (p24 = 9345 = <^i5 = 0, and 
^35 = ^62 = where ip = fcfasin^^. Since the lattice 
has the symmetry of point group Csi,, the vector potential Ai 
is rotated by ±|7r to obtain the other two vector potentials. 
Then the total accumulated phases along the nearest-neighbor 
directions are found to cancel out because of the symmetry 
of honeycomb lattice. However, the total accumulated phases 
for the next-nearest-neighbor hopping along the arrowed di- 
rections of the dashed lines in Fig. la are just Lp. Therefore, 
the total vector potential and magnetic field can be written as 

A = ?ifcf'[sin(2yfc^) +cos(V3.Tfc^)sin(yfc^)]e2: 



B 



— V^?ifcf sin('\/3xfc^ 
-2hk^k^ [cos(2yfc2 ) 



cos(yfc2)e.y (2) 
f 2 cos(%/3a;fc|') cos(?/fc^)]e^. 



The contours of the magnetic field are plotted in Fig. Id, in 
which the red lines indicate where the magnetic field van- 
ishes. The total magnetic flux through each hexagon vanishes, 
as the Peierls phase accumulated along its edges is zero. As 
a consequence, the total Hamiltonian of this cold atomic sys- 
tem can be described by Eq. (1). With Fourier transforma- 
k Ek and b, = ^ Ek ^'""^'h, the 



tion a, 



Hamiltonian of the system can be written by using "spinors" 

{ak,bkY as 



Hk = hoik) + hi{k)a^ + /i2(k)a2 + h3{k)a^ 



(3) 



where /io(k) = 2t' cos cos(k-ai), hi(k) =t^ -cos(k- 
s,), /i2(k) - iE.sin(k • s,), and haik) ^ M + 
2t sini^^,-sin(k • a^). (i = 1, 2, 3) are the three vec- 
tors pointing from a B site to its three nearest neighbors. The 



3 



vector /i(k) = [/ii(k), /i2(k), /i3(k)] is an effective magnetic 
field of the "spinors". The energy spectra are E = /io(k) ± 
|/i(k)| and the energy gap is {h^l, where ^13 = /i3(if+)( or 
h3{K^)) with K± = ±|2-(1,0). Ifhs = 0, the conduction 
band touches the valence band at two Dirac points A% and 
. The famous TKNN index or Chern number for this sys- 
tem is given by C = ^ J^^^ cPke^'^h ■ [dk^h x dk„h) 112611 
with /i(k) as the unit vector of h{k). It is demonstrated that 
the gauge invariance C can only take integer values and the 
quantum Hall conductivity is proportional to the Chern num- 
ber 10] ■ For the Haldane's model, different phases of the sys- 
tem can be characterized by different values of C. The phase 
diagram is depicted in Fig. (2a)||5l], in which the solid line is 
the critical boundary between the normal insulator with C = 
and the Chern insulator with C = ±1. It is notable that the 
phase (f ^ TT tan 9 can be controllable by simply choosing 
the laser angle 6, while the energy difference AI between sub- 
lattice can be tuned by the phase of laser beams. With such 
controllability, it is promising to realize the exotic topological 
phase transition between different Chern numbers in Fig. 2a. 

We now turn to establish a direct connection between the 
topological Cherm number and the atomic density, noting that 
the latter can be detected with density profile measurements 
typically used in atomic systems. We first develop a Green's 
function method to calculate the atomic density. The system is 
actually described by a Dirac-like Hamiltonian which can be 
obtained by expanding Eq.Q around two Dirac points K± = 
±||(1,0). By the substitution of k ^ K± + p, we have 
H± ~ ^vppia^ — + »Ti±(T3 at K±, respectively, 

where vp ~ is the Fermi velocity and m± = M ± 

sin (p. Under unitary transformation a^H-o^, we can 
write the Hamiltonian in a more symmetric form 



H± = -vppia^ - vpp2(7^ 



(4) 



where the notation m = ±m± is introduced for simplicity 
(i.e., we temporally omit subscripts ± below). 

It is observed that the conductivity axy (Chern number) of 
the system is related to the atomic density p according to the 
Streda formula axy = dp/dB\i^i^T once an additional uniform 
magnetic field B is applied. Such a magnetic field can be sim- 
ulated by rotating the optical lattice at a constant frequency 
uj = eB/2m. We choose the vector potential as ^0 = -^1 = 
and A2 = Bx. By using the substitution p ^ p + eA with 
e > 0, Eq.dHi can be solved in the real space. The eigenener- 
gies of the Hamiltonian can be obtained as jstl 

" \ ±^^771^ + 2nnvj,\eB\ n=l,2,3... ' ^ ' 

and the degeneracy of each LL is |/B|/0o per unit area with 
(f>Q the flux quantum. The density p in terms of the Green's 
function for the Dirac Hamiltonian is given by 

(0 + m)G = l, p=-Tr[7"G(x,x (6) 

where _0 ~ j'^Dr, with Dq = hdo — yU, Di = hvpdi, D2 — 
hvpd2 + ievpBx, 7° = a^, 7^ = a^, and 7^ = —a^ in 




FIG. 2: (a) Phase diagram of the system, where the phase boundary 
(solid line) corresponds to =0. The system behaves like a nor- 
mal insulator when C = and a Chern insulator when C = ±1. 
(b)-(e) Charge density in units S/(/)o for |m_| = 0.5|?n+| and 
eBhvp = 4|m+| as a function of normalized chemical potential 
fi/\m+ 1 (corresponding to a reseated atomic density profile in a trap) 
for four different cases: (b) m+ < < (c) m+ > > m_, 
(d) m+ < 771- < 0, and (e) ?n+ > m_ > 0. In (b) and (c) the shift 
of charge density at /i = is present so that the Hall conductivity is 
nonzero. In (d) and (e), the shift of the atomic density is absent and 
so the Hall conductivity vanishes. 



Euclidean space. From the standard Green function approach, 
the atomic density is expUcidy obtained as 



P± 



B 

I — |sgn(/i) <j int 



2nvj,\eB\ 



-}eM-\m±\ 



±- 



B 



m± 



200 \md 



e(|m±|-H), 



(7) 



where Q stands for the unit step function and int[a;] means the 
largest integer less than x. The second term of Eq. (|7]) is the 
atomic density induced into vacuum as /i ^ by the uniform 
magnetic field. It is of parity anomaly since its corresponding 
Hall cuiTent is independent of the magnetic field after multi- 
plying the drift velocity £ /B, where £ is the electric field. 

The total atomic density is given by the sum of the den- 
sities of the two components p = /?+ + /?_. At p, = 0, 

the Hall conductivity at B = can be obtained from the 

2 

density by using the Streda formula as a^y ~ Cj^, where 
C = ■^[sgn(m-|_) — sgn(m_)] is the Chern number To show 
how to detect the Chern number of the system, we consider 
a finite magnetic field B. The calculated density p in unit of 
\eB\/4>Q is plotted as a function of the normalized chemical 
potential ^/ 1 m+ 1 (for |m_| = 0.5|m+| wdehv'^B — 4|to+|) 
in Fig. 2. It is essential that the spatial density profile p{r) 
is uniquely determined by the function p{p/\m+\) in the lo- 
cal density approximation, which is typically well satisfied for 
trapped fermions. Figures 2b-2e stand for four different cases. 
The plateaus in the atomic density have one-to-one correspon- 
dence to the plateaus in the Hall conductivity due to the finite 
magnetic field S > 0. We here focus on /i = 0, which is of 
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our main interest. For m+ < < to_, which corresponds to 
C = —1, the atomic density p = -~B/4>o < 0, as shown in 
Fig. 2b. For m_|_ > > ?ti_, which corresponds to C = 1, 
the density p = B/<f>o > 0, as shown in Fig. 2c. For the other 
two cases to+ < to_ < and > to_ > 0, to+ and m_ 
have the same sign, corresponding to C = 0, and the density 
p = 0, as seen from Figs. 2d and 2e. Therefore, a simple di- 
rect relation between the Chern number and the equilibrium 
atomic density is established as 



C 



ho/B. 



(8) 



The important relation (O actually provides us a feasible 
way to experimentally detect the Chern number C in different 
phases. In the absence of B, the density of the cold atoms 
at /i = is first measured, which is denoted as po. Then 
the optical lattice is rotated to generate the effective uniform 
magnetic field B, and the new density of the cold atoms pi 
is measured. If pi > po, the system is in a Chern insulator 
phase with Chern number C = 1. If pi < po, the system is 
still a Chern insulator with C = —1. However, if pi — po, 
the system behaves like a normal insulator with Chern number 
C = 0. Since the density difference is actually quantized in 
units B/(j>Q, the above method could be rather robust. 

Finally, we briefly address an alternative approach to real- 
ize the Haldane's QHE. The fermions we discussed are in the 
s-band of the honeycomb lattice. As for fermions in the p- 
orbital bands, a Haldane's quantum Hall model without LLs 
can also be implemented by rotating each optical lattice site 
around its own center i27ll . Nevertheless, how to detect such 
p-band QHE is still desirably awaited. 

In summary, we have shown that the Haldane's QHE model 
can be realized by using ultracold atoms in an optical lattice. 
We have established a relationship between the Hall conduc- 
tivity and the equilibrium atomic density, which provides a 
feasible way to experimentally detect the Chern number C in 
different phases. 
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Note added - After this work was completed we became 



aware that a relation between the density profile and the 
Hall conductivity in conventional QHE was also addressed in 
Ref . ll28ll . We thank Dr. H. Zhai for bringing our attention to 
that work. 
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